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Abstract Time-of-Use (TOU) electricity pricing pro¬ 
vides a new opportunity for industrial user to reduce 
their power costs and improve the efficiency of power 
energy. Although many methods for Economic Load 
Dispatch (ELD) under TOU pricing in continuous pro¬ 
cess industry have been proposed, there is difficulty in 
batch-type process industry because of that the power 
load units are not fixed but closely related with pro¬ 
duction planning and scheduling. In this paper, for hot 
rolling, a typical batch-type and energy intensive pro¬ 
cess in steel industry, a multi-objective production schedul¬ 
ing optimization model for ELD is proposed under TOU 
pricing, in which the objective is to minimize the power 
costs on the premise of ensuring the product quality. A 
NSGA-II based production scheduling algorithm is pro¬ 
posed to generate Pareto-optimal solutions, and then 
the TOPSIS based multi-criteria decision-making is per¬ 
formed to recommend an optimal solution to facilitate 
filed operation. Experimental results on practical pro¬ 
duction data show that the proposed method cut down 
the power costs by creating load units corresponding 
to electricity price and shifting loads to avoid on-peak 
time periods. 
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1 Introduction 

Time-of-Use (TOU) electricity pricing, a practical de¬ 
mand response program that aimed to improve the peak 
load regulation ability of power grid, have been imple¬ 
ment by many power supplier now. It provide a huge op¬ 
portunity for electricity users to implement Economic 
Load Dispatch (ELD), which means cutting down their 
electricity costs by reducing the power loads at on-peak 
periods or shifting the loads from on-peak to off-peak 
periods. 

Optimizing power costs under TOU pricing can be 
quite different from optimizing energy consumption, the 
industrial users need to adjust their production schedul¬ 
ing according to time-dependent electricity tariffs to re¬ 
duce their power costs. In recent years, the methods 
that taking advantages of TOU pricing to implement 
ELD for industrial users have become a hot spot area. 
Shrouf et al. [1] proposed a single machine scheduling 
problem model, in which each time period has an asso¬ 
ciated price, the objective of model is to minimize power 
costs while considering traditional scheduling perfor¬ 
mance measures. Fang et al. P] also considered the job 
scheduling on a single machine to minimize the total 
power costs under TOU pricing and proposed the prob¬ 
lem solving algorithms for uniform-speed and speed- 
scalable machine environments respectively. Mitra et 
al. [3j formulated a mixed integer linear programming 
model for continuous process industry that allows op¬ 
timal production planning and provided an case study 
for the time horizon of one week and hourly chang¬ 
ing electricity prices, furthermore they improved the 
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Fig. 1 Process flow of hot rolling production 


model with integration of the operational and strate¬ 
gic decision-making [I]. Ashok [5] presented a theoret¬ 
ical model for batch-type load of process industry and 
proposed an integer programming method to resched¬ 
ule their operations to reduce the electricity costs un¬ 
der time varying electricity tariffs, but what model he 
presented is an abstract theoretical model and difficult 
to be applied to production directly. Wang et al. [B] 
proposed an optimization model with both the consid¬ 
eration of power generation scheduling and batch pro¬ 
duction load scheduling in the iron and steel plant to 
minimize the electricity cost, which is convinced to be 
effective under TOU tariff, but in which the production 
load units in their model are determined by fixed pro¬ 
duction planning and scheduling, the results can not 
always be optimal. 

From above analysis, we know that there are great 
potential for ELD under TOU pricing based on pro¬ 
duction scheduling. Most of the existing literatures are 
focused on some specific problem in continuous process 
industry or the abstract problem model, for ELD in 
batch-type process industry, there is difficulty that the 
power load units are not fixed but closely related to the 
results of production planning and scheduling. 

Hot rolling, a typical batch-type and energy inten¬ 
sive process in steel production with the characteristics 
of strong schedulability, has become an important as¬ 
pect on production organization and energy saving in 
steel industry [7]. In general, the basic process flow of 
hot rolling production can be illustrated as Fig. [T] Hot 
rolling is mainly organized and carried out by the batch 
scheduling program in steel mill, the primary task of the 
program is arranging and sequencing slabs into rolling 
units with the objective of smoothing the jumps among 
adjacent slabs in width, gauge, and hardness, which will 
directly affect the product quality. Hot rolling produc¬ 
tion scheduling has attracted attention from academics 
and industry areas for a long time. The earlier proposed 
method by Kosiba et al. thought the steel production 
scheduling as a discrete event sequencing problem, and 
formulated it as traveling salesman problem J5|. Lopez 


et al. [5] modeled the problem as a generalization of 
the prize collecting traveling salesman problem with 
multiple and conflicting objectives and constraints, and 
proposed a heuristic tabu search method to determine 
good approximate solutions. Tang and Wang [10. pro¬ 
posed a modified genetic algorithm to solved the prob¬ 
lem based on the multiple travelling salesman problem 
model. Chen et al. m formulated a nonlinear integer 
programming model and Kim [121 gave some correction 
to the model, furthermore, Alidaee and Wang [T3] pro¬ 
posed a corrected integer programming formulation and 
reduced the variables needed in the model. For most of 
the models are single objective or transformed models 
based on the weighted-sum approach, Jia et al. [Mj for¬ 
mulated the problem as a multi-objective vehicle rout¬ 
ing problem with double time windows model and pro¬ 
posed a decomposition-based hierarchical optimization 
algorithm to solve the model. Soon after, he proposed a 
P-MMAS algorithm to solve the problem, multi-criteria 
decision-making is performed to recommend the opti¬ 
mal solution from the Pareto-optimal solutions US- 

Due to high energy consumption and rising energy 
costs in hot rolling production m, energy saving has 
also been considered except the traditional objective 
mentioned above. As shown in Fig. [TJ the steel slabs 
need to be heated to high temperature before being 
rolled, the total energy consumption for heating is af¬ 
fected by the batch scheduling. Considering that high 
Direct Hot Charge Ratio (DHCR) in steel production 
has significant benefits such as energy cost reduction, 
some scholars made great efforts to improve the DHCR 
of slabs while performing batch scheduling [mim Be¬ 
sides that, optimization of rolling schedule by adjusting 
reduction ratio of slabs between the rolling passes, an¬ 
other way to reduce the power consumption that used 
to drive the rolling motor, is also proposed msm 
and has taken a certain effect. 

As mentioned above, almost all the existed meth¬ 
ods for hot rolling production scheduling are concen¬ 
trated on internal production organization. Although 
some technical means have been proposed and applied 
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Fig. 2 diagrammatic sketch of batch scheduling: (a) rolling unit, (b) rolling batch. 


to achieve energy saving in some studies, the potential 
would have been exhausted because of the restriction of 
equipment and technology constraints. In this context, 
the methods that utilizing the favorable external en¬ 
vironment should be explored for energy saving. TOU 
pricing provides an opportunity to reduce their power 
costs, but until now the papers to implement ELD un¬ 
der TOU pricing for hot rolling production are few pub¬ 
lished. 

This paper consider the Hot Rolling Production Sche¬ 
duling Problem (HRPSP) as a mixture of batch schedul¬ 
ing problem and time-dependant job-shop scheduling 
problem, the rolling units, which are seem as power 
load units in this paper, are planed and scheduled ac¬ 
cording to the TOU electricity prices. The primary ob¬ 
jective of this model is to minimize the power costs 
while considering the traditional objective of minimiz¬ 
ing penalties that caused by jumps between adjacent 
slabs. A multi-objective optimization model and algo¬ 
rithm are proposed to solve the problem, and a multi¬ 
criteria decision-making is performed to recommend an 
optimal solution to facilitate field operation. 


2 Problem description and formulation 

HRPSP is an extremely complex problem that has been 
convinced to have significant influence on product qual¬ 
ity, production efficiency and energy consumption. In 
this paper, we consider the HRPSP as Hot Rolling Batch 
Scheduling Problem (HRBSP) combined with Job-shop 
Scheduling Problem (JSP), in which the HRBSP is fo¬ 
cused on how the rolling units be scheduled and the 
JSP is concentrated on how and when the rolling units 
be produced. 


2.1 Problem description 

Hot rolling batch scheduling is a key process in hot 
rolling production. In this paper, we follow the scenario 
as shown in Fig. [2] to describe it as follows: The task 
of HRBSP is selecting, grouping and sequencing can¬ 
didate slabs into rolling units under the constraints of 
production capacity and rolling rules. As shown in Fig. 
[2l a rolling unit that has a coffin-shaped width profile 
consists of a warming-up section and a coming-down 
section. In the previous one the slabs are arranged from 
narrow to wide to warm up the rolls, and in the later one 
the slabs are scheduled from wide to narrow to avoid 
marks on the surface of the coils. The major part of a 
rolling unit is the coming-down section, and the quality 
of rolling unit mainly depends on the sequence of slabs 
in it. In most cases, the warming-up section is just a 
minor part that can be determined manually. 

There are several constraints that restrict the rolling 
batch scheduling, in which the most important one is 
smoothing jumps between adjacent slabs in width, gauge, 
and hardness. Other constraints, such as the cumula¬ 
tive rolling length of slabs in a rolling unit, the contin¬ 
uously rolled length of slabs with the same width and 
so on, should also be considered to ensure the prod¬ 
uct quality and production capability. For hot rolling is 
an key energy intensive process in steel industry, many 
approaches, such as optimization of batch scheduling 
with the objective of improving hot charging ratio and 
optimization of reduction schedule to reduce power en¬ 
ergy for driving the rolling motor, have been proposed 
to achieve energy saving. In smart grid, TOU electric¬ 
ity pricing, which have been one of the most common 
implementation of demand response programs [22|, pro¬ 
vide a new opportunity for the steel mill to achieve ELD 
in hot rolling production, which means cutting costs by 
shifting loads among different time periods. 

As shown in Fig. [3j a whole day is split into four 
types of time periods: on-peak, mid-peak, flat-peak and 
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Fig. 3 Relationship between production scheduling and power costs under TOU pricing 


off-peak periods, different types of time periods are as¬ 
sociated with different electricity prices. We can also see 
that the power costs for each rolling unit, which are not 
only determined by the quantity of the power demand 
for each slab but also depended on the corresponding 
TOU electricity pricing, should be accumulated piece- 
wise during the processing time. 

Compared to the flat electricity pricing, the objec¬ 
tive of ELD is to minimize the total power costs, which 
consists of the charges for power consumed and the ad¬ 
ditional operating costs due to shifting of loads, if any. 
In this paper, we assume that the rolling units can be 
scheduled freely, and therere no additional operating 
costs due to shifting of loads. In this context, the rolling 
production arranged in off-peak periods is encouraged 
and that in on-peak periods is deprecated to implement 
the peak load shifting. At the same time, the idle time 
should be allocated in the periods with high price as 
far as possible. In addition, we should know that the 
scheduling on fixed jobs are not always optimal, so the 
scheduled jobs, which means the rolling units generated 
by hot rolling batch scheduling, should be created with 
association to their operation time. 

Finally, the problem is turned into optimal schedul¬ 
ing for minimizing the power costs that determined by 
batch scheduling solution and the job-shop scheduling 
solution under specified electricity pricing, while the 
traditional objective that smoothing the changes be¬ 
tween adjacent slabs should not be ignored to ensuring 
the product quality. 


2.2 mathematical formulation 

As in literature [15], the basic model of the HRBSP can 
be interpreted as a Price Collecting Vehicle Routing 
Problem (PCVRP), which is an extension of VRP that 
well-known as an NP-liard problem. In the PCVRP 
model, it can be considered that each rolling unit is 


a vehicle within limited capacity and each slab is a cus¬ 
tomer that can be visited at most once. Suppose that 
there are n slabs to be scheduled into m rolling units, 
the objective of the problem is to determine m routes 
(rolling units) to minimize the total distance traveled 
(penalties caused by jumps between adjacent slabs). 

Considering the relationship between the slab pro¬ 
cessing sequence and the power costs under TOU pric¬ 
ing as shown in Fig. [HI the HRPSP can be formulated as 
a variant of PCVRP, in which the notations are defined 
as follows: 

T: the sequence of time periods, T = {1,2,..., t}, 
where t is the number of time periods, a corresponding 
electricity price is associated to each period; 

A i'. the duration of the ith time period, where i £ T; 

7 r,: the electricity price at the ith time period, where 
i £ T; 

N: the sequence of all candidate slabs in the slab 
pool, N = {1,2,..., ?7.}, where n denotes the quantity 
of slabs; 

M: the sequence of rolling units, M = {1, 2,..., m}, 
where variable m is the quantity of rolling units; 

Pij-. the penalties to roll the jth slab immediately 
after the ith slab, where = pfj + pf- + Pij, and 
the variables pfj , pC and p^ respectively represent the 
penalties that caused by jumps from slab i moving to 
slab j in width, gauge, and hardness; 

If the length of the ith slab, where i £ N; 

wf the width of the ith slab, where i £ TV; 

L: the lower bound of the cumulative length of slabs 
that scheduled in a single rolling unit; 

U : the upper bound of the cumulative length of slabs 
that scheduled in a single rolling unit; 

R : the upper bound of the length of continuously 
rolled slabs with the same width in a single rolling unit; 

Wf the power consumption of the ith slab during 
rolling process, where i £ N; 

Pf the processing time needed for the ith slab, where 

i £ iV; 
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vp. the idle time allocated to the jth rolling unit 
before the production, where j G M; 

TS: the total time that can be allocated for produc¬ 
tion. 

To identify how the slabs are scheduled in rolling 
units and when the slabs are processed, four binary 
decision variables are defined as 

yk = I 1 if * h 

I 0 otherwise. 


subject to 

Y4 = vt,jeN,keM 

i£N 


Y x io = y1 ?' 

jeiv 


where i b j means the jth slab is scheduled in the jth 
rolling unit; 


fceM 


x k = 


1 if Vi = 1, Vj = 1 and j y i, 
0 otherwise. 


where j y i means the jth slab is rolled immediately 
after the jth slab; 


s k . 


\ if x k j = 1 and Wi = Wj, 
0 otherwise. 


Y 4 ■ Vj -lj<R,ie N, k G M 

j£N 

L < Y Vi ■ U < u > k 6 M 

iGN 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 


which means the jth slab is rolled immediately after the 
jth slab in the fcth rolling unit with the same width; 

f if "YhaKj — Pi — Yla<j 
0 otherwise. 

where Y^ineqi ^ a means the duration of time periods 
subjected to the inequality Ineql , p[ means the pro¬ 
cessing start time for the jth slab, which is not only 
determined by decision variables xj and yj , but also de¬ 
pended on the processing time of rolling units, the idle 
time allocated for rolling units and duration of TOU 
time periods. Furthermore, p[ can be expressed as 

Pi = Y Y y “' Pa + Y Va+ J2 X if>(k) ■ Pen 

k<5 a£N a<5 k<6(i) a£N 

where the variable S represents the sequence number 
of rolling unit that the jth slab is scheduled, 0(i) is a 
simple index function to return the position of the jth 
slab in the rolling unit, Y^ae N Pa 'Pa is the cumulative 
processing time before the <5th rolling unit, Y^ a <s Va 
the cumulative idle time allocated before the jth rolling 
unit, J2k<e(i) T,acN X ie(k) 'Pa represents the cumula¬ 
tive processing time before the jth slab in the <5th rolling 
unit. 

Based on the above notations and decision variables, 
a multi-objective hot rolling production optimization 
model can be formulated as 

min /i = Y l 71 * ■ Y12 Wi ' d i J w 

fceT \ ieJVjeT J 

min h=YYY ( P b • x d k ) (2) 

fe£M i€N j£N 


YdTs-YYy*-^ ( 8 ) 

i£M fceM i£N 

where the objective CD means to minimize the total 
power costs in hot rolling production to implement eco¬ 
nomic load dispatch, and the objective © is the tra¬ 
ditional objective to ensure the product quality, which 
means to minimize the total penalties caused by jumps 
between adjacent slabs. Constraints © and © spec¬ 
ify the sequence of slabs in a rolling unit. Constraint 
© ensures that each slab can be scheduled only once. 
Constraints © restrict the cumulative length of contin¬ 
uously rolled slabs with the same width in each rolling 
unit. Constraint © is the constraint of rolling mill pro¬ 
duction capacity, which restricts the lower and upper 
bounds of cumulative length of slabs in each rolling 
unit. Constraint © means that the total idle time allo¬ 
cated for rolling units cannot be greater than the mar¬ 
gin of production capability. 


3 Multi-objective production scheduling 
method 

In this paper, the production scheduling method con¬ 
sists of two stages to implement ELD. In the first stage, 
objectives shown in Eq. © are optimized simultane¬ 
ously, and then a set of Pareto-optimal solutions is gen¬ 
erated by the multi-objective optimization algorithm. 
In the second stage, the TOPSIS based multi-criteria 
decision-making is performed to recommend an opti¬ 
mal solution to facilitate field operation. 
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3.1 Multi-objective production scheduling algorithm 

NSGA-II is a classical multi-objective evolutionary op¬ 
timization algorithm that proposed in 2000 by Deb et 
al. [23]. Because of the advantages of low computa¬ 
tional expenditure, elitist strategy and less parameter 
sharing approach, NSGA-II has been widely used to 
solve the combinatorial optimization problems in en¬ 
gineering, such as the hydro-thermal power scheduling 
problem m , job sequencing problem m and flow-shop 
scheduling problem 20 . 

In this paper, a NSGA-II based Multi-objective Pro¬ 
duction Scheduling Algorithm (MOPSA) is developed 
to solve the HRPSP, and the overall procedure of MOPSA 
is inherited from NSGA-II. Some personalized changes 
are made to instantiate the algorithm, in which the 
most import things are designing customized chromo¬ 
some code and genetic operators to adapt specific prob¬ 
lem. 

For containing information about the batch schedul¬ 
ing and job-shop scheduling, a hybrid chromosome code 
consists of two parts as shown in Fig. [4] is designed. 
The first part is a natural number sequence C that can 
be transformed to a two-dimensional matrix B through 
the code mapping procedure, the matrix B represents a 
batch scheduling solution and the element bj in B is the 
original sequence of the jth slab in the «th rolling unit. 
For each i. if the minimal of j is found while bij = 0, 
it can be resolved that the last slab in the ith rolling 
unit is bij- 1 . The second part is a floating number se¬ 
quence V representing the idle time allocated for job- 
shop scheduling, where job means the production of 
rolling units. 



Fig. 4 Relationship between production scheduling and 
power costs under TOU pricing 


According to above description, the hybrid chromo¬ 
some code G can be formulated as 

(G = (C,V) 

< c = (ci,c 2 ,...,c m x n) , 

{V = (U1,U 2 , . . . ,Vm) 

where the element Cj in C is a natural number that 
ranged from 1 to to x n, to is the quantity of rolling 
units and n is the quantity of slabs to be scheduled, 
any two number Cj and Cj are associated with different 
values, Vi in V represents the idle time allocated to the 
ith rolling unit before rolling production. 

Detailed steps of the code mapping procedure as 
mentioned above are listed as follows: 


Step 1: Set fi{i = 1, 2,..., n) to 0, where is a flag 
and fi = 1 represent the itli slab has been arranged into 
rolling units; for the fcth (k = 1 , 2 ,..., to) rolling unit, 
set nurrik = 0 , which means that the quantity of slabs 
in it is 0 , set dk = 0 , which means that the accumulative 
rolling length of it is 0 , set qk = 0 , which means that 
the continuously rolled length of slabs with same width 
is 0 ; set a loop variable j = 1 ; 


Step 2: Confirm the variables s and k that correspond¬ 
ing to the natural number Cj , by which it can be deter¬ 
mined that the sth slab are scheduled in the kth rolling 
unit, where s and k can be calculated by 



x TO 


and 



Step 3: Check if the condition f s = 0 is satisfied: 

i. If satisfied, it means that the sth slab is an unsched¬ 
uled slab. Then if w s ^ w' k , set qk = 0, where w s is 
the width of the sth slab and w' k is the width of the 
latest appended slab in the fcth rolling unit. Further¬ 
more, if dk + l s < L and qk + l s < R, put the sth 
slab into the kth rolling unit, set nurrik = numk + 1 , 
dk = d k + l s , qk = qk + l s and f s = 1 , the matrix 
B(= [bj]) that represents slabs in rolling units will be 
updated as b k ^ numk = s; 

ii. Otherwise, go to step 4; 


Step 4-' Update variable j = j + 1, go to step 2 to repeat 
the above operations until j=mxn+l; 
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Step 5: Check if all the /* = 0(i = 1,2 .,n) are sat¬ 
isfied: 

i. If satisfied, it is means that all the slabs are all allo¬ 
cated to rolling units with the constraints, and the ma¬ 
trix £?(= [bij]) represents a feasible solution of rolling 
batch; ii. Otherwise, a large number should be assigned 
to the fitness function /i to avoid the chromosome be 
selected into the new population in the next selection 
operator. 


We know that the constrained optimization is usu¬ 
ally more difficult than the unconstrained optimization. 
The benefits of hybrid encoding and mapping procedure 
in this paper are not only containing complete informa¬ 
tion of the production scheduling but also handling the 
constraints. From step 5, it can be seen that all the con¬ 
straints from Eq. © to © in Sec. IQ are satisfied in 
the finally accepted feasible solution, which is helpful 
to reduce the difficulty for problem solving. 

Then in order to instantiate the MOPSA algorithm, 
customized genetic operators should be defined to match 
the hybrid chromosome code. For genetic algorithm, the 
most important genetic operators are selection, crossover, 
and mutation, which are briefly introduced as below. 
Selection operator, selecting individuals from popula¬ 
tion, is done based on the front rank of individuals by 
non-dominated sorting. For individuals with the same 
rank, which have the larger crowded distances will be 
adopted preferentially. In order to generate offspring 
from parent chromosomes, Partially Mapped Crossover 
(PMX) [27] and Scramble Sub-list Mutation (SSM) [28] 
are adopted to operate on the part of C in the chromo¬ 
some code. The PMX operator is performed on two par¬ 
ent chromosomes: randomly select two crossover points 
k\ and , swap the gene codes in the range [k\ , /C2] , af¬ 
ter that, replace the other gene codes out of the range 
[ki , ^2] according to the mapping relationship that de¬ 
termined by the middle section. Unlike the crossover 
operator, the SSM operator is performed on single par¬ 
ent chromosome: randomly select two positions p\ and 
P 2 that separated less than a fixed length in the chro¬ 
mosome code, and then rearrange the gene codes be¬ 
tween [pi, P 2 ] • After crossover or mutation operators, 
the update on the part of V in chromosome code, which 
means idle time allocation for rolling units, should be 
performed immediately. 

The detailed steps of idle time allocation as men¬ 
tioned above for rolling units based on the part of C in 
the chromosome code are listed as follows: 


Step 1: Confirm the electricity price 7 r® while the pro¬ 
duction started and price 7 r? while the production com¬ 
pleted for the ith rolling unit, where i £ M; 

Step 2: Create a random floating number sequence V 
that represents the idle time allocated for rolling units 
before production. For elements V, the constraint as 
Eq. © in Sect. 12.21 must be satisfied; 

Step 3: Sort the time periods in descending order based 
on the electricity price, after that a new set of time 
periods T' = ■. ■ ,t' t ) is generated, in which the 

price associated with t' k is 7r(,; define a loop variable j 
and set j = 1 ; 

Step 4 : Adjust the idle time allocation for rolling units. 
Suppose that the ith rolling unit is started from time 
period i', if irf < 7 r' and Vi +1 > 0 , then set Vi -|_i = 
0 , Vi = Vi + |_i; suppose that the ith rolling unit is 

completed in time period f', if 7r| < 7 r' and Vi > 0, then 
set Vi = 0, v i+ i = v i+ i + Vi ; 


Step 5: Update variable j = j + 1, go to step 4 to 
repeat the operation for the left time periods until j = t , 
which represent the adjustment of idle time allocation 
is completed. 


3.2 TOPSIS based multi-criteria decision-making 

As MOPSA generate more than one Pareto-optimal so¬ 
lution, in order to facilitate field operation, only a few 
solutions should be accepted. In this paper, Technique 
for Order Preference by Similarity to an Ideal Solution 
(TOPSIS) :25], a widely used multi-criteria decision¬ 
making method to identify solutions from finite alter¬ 
natives, is adopted as the method to select a recom¬ 
mended optimal solution. 

Detailed steps of the TOPSIS based multi-criteria 
decision-making for HRPSP are listed as follows: 

Step 1: The decision matrix X can be expressed as 

X12 ’ 

X22 

. 5 

_ ^ml Xm2 J 

where X is a two dimensional matrix with the size of 
mxn, which means that therere m solutions generated 
by the multi-objective algorithm and two objectives for 
the HRPSP, element Xij in A' is the value of the jth 
objective with respect to the ith solution. Then the 


X = 


£11 

X21 
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normalized decision matrix Z(= [. Zij ]) can be calculated 
according to 



Step 2: Multiply the normalized decision matrix by its 
associated weights to calculate the weighted normalized 
decision matrix V(= [vij}), in which Vjj is calculated as 

Vij = wj ■ z ih 

where Wj is the weight associated with the jth objec¬ 
tive, and in this paper, w\ and w 2 are set to 0.4 and 0.6 
respectively with the preference of the two objectives. 


Step 3: Identify the the ideal solution s + and the nadir 
solution s~ of each objective according to the following 
equations: 

s + = (s^,4) , 

{ max Vij if fj is benefit-oriented, 

l<i<m 

min Vij if fj is cost-oriented. 

1 <i<m 

S = ( S 1 ) S 2 ) J 

! min if fj is benefit-oriented, 

1 <i<m 

max v^ if fj is cost-oriented. 

l<i<m 

It should be known that both of the objectives in 
HRPSP are cost-oriented, which is said to find the min¬ 
imum of objective functions. 


Step 4 : Measure the distance df and d i of the ith so¬ 
lution from the ideal solution s + and the nadir solution 
s" by 


df = 

d~ = 


2 

a X] ( Vi i 

\ i=i 


2 

a ( u b 
\ J= 1 


- s t) 2 ’ i= X ’ 2 ’'' ■ > m > 


_ v 2 

- Sj ) ,i=l,2,...,m. 


Step 5: Calculate C* that represents the relative close¬ 
ness of Ah solution with respect to the ideal solution 
according to 


G* = 


(d~ + d+) 




After completing the above steps, the decision-making 
can be finally performed on the Pareto-optinral solu¬ 
tions according to the sequence that determined by 
C*{i = 1,2,...,to) in descending order, which solu¬ 
tion that owns the maximal relative closeness will be 
selected as the recommended optimal solution. 


4 Experimental procedure and results 

To evaluate the effectiveness and advantages of the pro¬ 
posed method, four groups of production data with dif¬ 
ferent slab quantities and processing time are collected 
from Xiangtan Steel Mill to be used in the experiment. 
As shown in Table[lJ each group of data represents a dif¬ 
ferent characteristic. For the production data, many va¬ 
rieties means that there are numerous variety of width, 
gauge and hardness for the slabs, and the penalty be¬ 
tween adjacent slabs will be larger, at the same time, 
full load means that the slab processing will take long 
time and the idle time to allocate for rolling units will 
be fewer. 

The data in Table [2] is the TOU electricity tariffs 
that actually performed in the steel mill. According to 
daily power load distribution, the whole day is split into 
eight periods that covered four types of time periods, 
and for each type of time period, a corresponding elec¬ 
tricity price is assigned. 

According to constraints of equipment and produc¬ 
tion capability, the lower and upper bound of the cumu¬ 
lative length of slabs that scheduled in a single rolling 
unit are respectively set to 5 and 10 kilometers, and 
the upper bound of the continuously rolled length of 
slabs with same width is set to 1 kilometers. For a spe¬ 
cific slab, the rolling length, processing time and power 
consumption can be obtained by the hot rolling process 
control system in steel mill. The penalties that caused 
by jumps between adjacent slabs in width, gauge and 
hardness are adopted by referencing to [5]. 

To obtain excellent performance and preferable re¬ 
sults, choosing suitable parameters values for the MOPSA 


Table 1 Production data description 


Group Slab 
Id quantity 


Rolling 

units 

quantity 


Characteristics 


1 

442 

8 

Many varieties, full load 

2 

418 

8 

Many varieties, not full load 

3 

441 

8 

Few varieties, full load 

4 

417 

8 

Few varieties, not full load 


Table 2 Time-of-use electricity tariffs 


Time 

period 

Time frame 

Electricity price 
/{CNY ■ kWh- 1 ) 

on-peak 

18:00-21:00 

0.878 

mid-peak 

08:00-11:00, 15:00-18:00 

0.778 

flat-peak 

07:00-08:00,11:00-15:00, 

1:00-22:00 

0.628 

off-peak 

00:00-07:00, 22:00-24:00 

0.428 
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Penalty 


Fig. 5 Pareto fronts for different groups of production data: 

algorithm is very important. Using parameter sensitiv¬ 
ity analysis based on a lot of test, a group of empir¬ 
ical values is determined, in which the probability of 
crossover and mutation are set to 0.4 and 0.6 respec¬ 
tively, the size of population is set to 50, and the max¬ 
imum iterations of algorithm is set to 2000. 

The proposed method is implemented in MATLAB 
7.0, and the experimental results are obtained on a PC 
with Pentium(R) 2.13 GHz dual core processor and 2 
GB memory. 


4.1 Experimental results 

To measure the effectiveness of the proposed method, 
the graphical approximate Pareto fronts returned by 
MOPSA are illustrated in Fig. [5l in which the sub dia¬ 
grams (a) to (d) are related to four groups of production 
data respectively. In each sub diagram, the horizontal 
axis denotes the total penalties caused by jumps and 
the vertical axis denotes the power costs in production. 
Because both of the objectives in our problem are cost- 
oriented, we know that the smaller of objectives, the 
better performance of the algorithm will be. 



Penalty 


(a) group 1, (b) group 2, (c) group 3, (d) group 4. 

As shown in Fig. [5] we can see that for each group 
of production data, a certain number of Pareto-optimal 
solutions are generated and uniformly scattered on the 
Pareto fronts, which means that the MOPSA has strong 
ability for searching optimal solution. For sub diagrams 
(a) and (b), the penalty scores are obvious larger than 
(c) and (d), that is because the varieties of slab in group 
1 and 2 are varied and abundant. At the same time, the 
power costs in (a) and (c) are larger than (b) and (d), 
which is because that the production load in group 1 
and 3 are fuller. Nevertheless, it should be noted that 
the true Pareto-optimal front is not known for the prac¬ 
tical production data, so the fronts in Fig. [5] are all 
pseudo-optimal Pareto fronts constructed by combin¬ 
ing all the non-dominated solutions and removing all 
the dominated solutions. 

In order to facilitate field operation, the TOPSIS 
based multi-criteria decision-making is performed to se¬ 
lect the only one solution from the Pareto-optimal so¬ 
lutions. As shown in Fig. [5] the dimension of two ob¬ 
jectives are different, which means /i is much bigger 
than y *2 and fi fluctuate in a smaller range than / 2 , 
which will affect the practical importance of objectives 
in TOPSIS sorting. In order to put the values of two 
objectives into similar scale, the transformation for el- 
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Table 3 Partial TOPSIS sorting results 


Solution Id 

Relative closeness 

fi (/CNY) 

h 

SI 

0.7772 

300600.61 

7097 

S2 

0.7758 

300414.79 

7424 

S3 

0.7458 

300079.79 

8081 

S4 

0.7418 

300296.46 

7949 

S5 

0.6845 

301898.17 

6654 

S6 

0.5627 

299376.89 

9916 

S7 

0.4514 

299368.24 

10899 

S8 

0.3155 

298841.41 

13330 


Table 4 

Comparison of scheduling solutions 


Group 

Proposed method 

Conventional method 

Id 

Power 

cost/(CNY) 

Penalty 

Power 

cost/(CNY) 

Penalty 

1 

300600.06 

7097 

309389.81 

6737.00 

2 

278606.13 

7049 

297575.21 

6265.00 

3 

302401.27 

3039 

310091.60 

2765.00 

4 

278330.86 

3300 

299708.13 

2923.00 


ements in decision matrix is performed as 

Xij = Xij — min Xij, 

1 <i<m 

where means the value of the jth objective with 
respect to the ith solution, and variable m is quantity 
of Pareto-optimal solutions. 

After the transformation, TOPSIS sorting is per¬ 
formed and partial results for the first group of pro¬ 
duction data are given in Tabic [3] According to the de¬ 
scending order based on relative closeness, the solution 
SI, which owns the maximal relative closeness 0.7772, 
is finally selected as the recommended optimal solution. 

Comparisons of scheduling results by two methods 
are given in Table [I] The method in m is adopted 
as the conventional method, in which the penalty score 
is minimized by GA algorithm and the power costs is 
never taken into consideration. It can be seen clearly 
that the penalty by the proposed method increased 
than that by traditional method but the magnitude is 
very small for any group of data. Compared to reference 
value in turn the penalty is still accepted in engineer¬ 
ing, and the product quality would not be affected obvi¬ 
ously, but at the same time, much lower power costs are 
obtained by the proposed method with the same power 
consumption, which indicating that ELD, the primary 
objective of the proposed method in this paper, is im¬ 
plemented by utilizing the TOU pricing in electricity 
market environment. 

4.2 Analysis of scheduling results 

To further confirm the effectiveness of the proposed 
method, the job scheduling results on the first group 


of production data is chosen as a scenario to perform a 
more detailed analysis. Values of rolling parameters by 
the two method are both given in Table [H from which 
we can see that the quantities of slabs in different rolling 
units are basically uniform, and the rolling length for 
each rolling unit is restricted to less than 10 kilometers. 
Besides that, some other parameters for each rolling 
unit, such as processing time, power demand, average 
power load, processing start time, processing complete 
time, and allocated idle time, are also given to describe 
the results of job scheduling, all of the above param¬ 
eters represent that the proposed method is effective 
while owning the better performance on objectives /i 
as shown in Table 0 ] 

Considering rolling units as production jobs and il¬ 
lustrating them in a double longitudinal axis diagram 
as shown in Fig. [G] in which the sub diagrams (a) - (d) 
are corresponding to group 1 to 4 of production data 
respectively. The horizontal axis in sub diagram repre¬ 
sents the rolling processing time, and the left and right 
longitudinal axes indicate the power load and electric¬ 
ity price respectively. The jobs scheduled by proposed 
method are represented by solid lines that started and 
ended with a circle, and the jobs scheduled by conven¬ 
tional method are illustrated by solid lines that started 
and ended with a square. From the figure we can iden¬ 
tify the production load and the electricity price at any 
time. 

As it can be seen, in any sub diagram, the heavy 
loads are created and mainly arranged in off-peak and 
flat-peak periods by the proposed method, while the 
light loads are arranged at on-peak or mid-peak periods. 
In addition, the idle time is allocated for jobs under 
constraints of production capability by the proposed 
method. For scenario in this paper, we can see that the 
idle time is allocated at 18:00 to 21:00, which is the 
obvious on-peak time period. 

Form above analysis, we can know that the advan¬ 
tages of the proposed method can be attributed to two 
aspects, one is the TOU pricing based power load cre¬ 
ating, which means the rolling units are not fixed but 
depended on their processing time and the correspond¬ 
ing electricity price, another one is load shifting, which 
means arranging the processing sequence and time for 
rolling units to avoid on-peak time periods. 

4.3 Discussion of peak load regulation 

From perspective of power grid, the intention of im¬ 
plementing TOU electricity pricing is encouraging the 
power users to participate peak load regulation, which 
means reducing power loads in on-peak periods or shift¬ 
ing power loads from on-peak to off-peak periods, and 
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Table 5 Detailed comparison of the scheduling solution for the 1st group of production data 


Scheduling 

method 

Rolling 
unit No. 

Slab 

quantity 

Rolling 

length 

/(km) 

Processing 

time 

/(h) 

Power 

demand 

/(MW-h) 

Average 
power load 
/(MW) 

Processing 

start 

time 

Processing 

complete 

time 

Allocated 

idle 

time 


1 

59 

9.12 

2.67 

59.71 

22.36 

00:00 

02:40 

0 


2 

50 

8.60 

2.67 

57.89 

21.68 

02:40 

05:20 

0 


3 

51 

9.97 

2.84 

60.96 

21.47 

05:20 

08:10 

0 

Proposed 

4 

59 

9.48 

3.01 

60.57 

20.12 

08:10 

11:11 

0 

method 

5 

57 

9.97 

3.11 

65.64 

21.11 

11:11 

14:18 

0 


6 

57 

8.71 

2.89 

58.37 

20.20 

14:18 

17:11 

0 


7 

50 

9.98 

3.15 

63.88 

20.28 

17:11 

20:20 

0 


8 

59 

9.99 

2.99 

64.72 

21.65 

21:00 

24:00 

0.66 


1 

46 

7.35 

2.47 

50.63 

20.50 

00:00 

02:28 

0 


2 

58 

9.94 

3.03 

65.08 

21.48 

02:28 

05:30 

0 


3 

57 

9.88 

2.95 

62.69 

21.25 

05:30 

08:27 

0 

Conventional 

4 

57 

9.98 

3.02 

63.27 

20.95 

08:27 

11:28 

0 

method 

5 

51 

9.05 

2.74 

57.83 

21.11 

11:28 

14:12 

0 


6 

60 

9.73 

3.16 

65.98 

20.88 

14:12 

17:22 

0 


7 

56 

9.93 

2.92 

61.93 

21.21 

17:22 

20:17 

0 


8 

57 

9.94 

3.05 

64.33 

21.09 

20:17 

23:20 

0 
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Fig. 6 Scheduling results on different groups of production data: (a) group 1, (b) group 2, (c) group 3, (d) group 4. 


is good for promoting the stability of power grid. In 
this section, we try to discuss the effectiveness of the 
proposed method on peak load regulation. 

The power load distribution among time periods by 
the two scheduling methods are provided in Fig. [TJ in 
which the comparisons for all groups of production data 
are illustrated by the sub diagrams (a) - (d) respec¬ 


tively. As it can be seen, for any group of production 
data, the average power load is greatly reduced in the 
on-peak periods and substantially increased in the last 
off-peak period, at the same time, the power loads in the 
first off-peak period increased slightly. Especially, for 
group 2 and group 4, in which the production data are 
charactered by not full production load, the power load 
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Fig. 7 Power load distribution among time periods for different groups of production data: (a) group 1, (b) group 2, (c) group 
3, (d) group 4. 


increased in the on-peak period and reduced in the last 
off-peak period significantly. Combined with Fig. 6, we 
can know that the above results are mainly attributed 
to power load shifting from on-peak to off-peak peri¬ 
ods. From the results in Fig. 7 and the above analysis, 
it can be concluded preliminarily that the ELD in this 
paper promoted the peak load regulation ability of the 
power grid, and the effectiveness of peak load regulation 
is majorly affected by the full condition of production 
load. The idler of production load, the more potential 
of the peak load regulation ability will be. 

Nevertheless, it should be noted that the ELD dont 
always promote the peak load regulation. Considering 
that a new TOU tariff, in which an on-peak period is 
closed to the off-peak period, is provides to users as 
shown in Fig. [51 and three jobs need to be completed. 
In order to minimize the power costs, the Jobl and 
Job3 are arranged as close as possible to the off-peak 
periods. In order to implement peak load regulation, the 
Job2 should be arranged in the mid-peak or flat-peak 
periods, but considering minimizing the power costs, 
the left Job2 would be arranged across the off-peak and 
on-peak period, in which most of the job processing 
time will be covered by the off-peak period. In this case, 



Fig. 8 Scheduling results under a new TOU tariffs 

the power load in the first on-peak period will not fall 
but rise, which deviate from the original intention of 
TOU pricing. This is a special case that arising from 
the lack of interaction between the power supply side 
and the demand side in electricity market, and it needs 
to be solved with the participation and collaboration 
of both sides, which will not be discussed in depth in 
this paper. In our method, we just perform the optimal 
production scheduling according to the provided price 
signals to reduce the power costs. 
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5 Conclusions 

This paper presented the challenge of energy saving in 
hot rolling production and formulated a multi-objective 
optimization model of HRPSP under TOU electricity 
pricing, in which the hot rolling batch scheduling and 
the time-dependant job-shop scheduling were both taken 
into consideration. Objective of the proposed model is 
to minimize the power costs while considering tradi¬ 
tional objective for measuring the scheduling perfor¬ 
mance. A NSGA-II based multi-objective production 
scheduling algorithm was proposed to solve the prob¬ 
lem, in which the chromosome code and genetic op¬ 
erators are custom-designed to match the HRPSP, and 
moreover, a TOPSIS based multi-criteria decision-making 
was performed to recommend an optimal solution from 
the Pareto-optimal solutions to facilitate filed opera¬ 
tion. Experimental procedure on practical production 
data was given out, the results and analysis showed 
that the proposed method can reduce power costs in 
hot rolling production on the premise of ensuring the 
product quality. Considering multiple production lines 
existed in most steel mills, the further study on HRPSP 
integrated parallel machine job-shop scheduling will be 
the next work, which is expected to be more difficult 
but have greater benefits. 
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